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ABSTRACT 

We present global radiative transfer models for heated relativistic jets. The simulations include all 
relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation 
to and beyond the photosphere of the jet. The transfer models are compared with three gamma- 
ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different 
spectra. The models provide good fits to the observed spectra in all three cases, and we obtain 
estimates for the jet magnetization parameter eb and the Lorentz factor T. In the small sample of 
three bursts, Eb varies between 0.01 and 0.05, and T varies between 400 and 1200. 

Subject headings: gamma-ray burst: general - plasmas - radiation mechanisms: non-thermal - radia¬ 
tive transfer - scattering 


1. INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) have been ob¬ 
served for more than four decades, and thousands of 
bursts have rather well measured spectra. No physical 
model has been demonstrated to systematically fit the 
observed spectra. Instead, the data are usual ly described 
using a phenomenological Band function (iBand et al.l 
Il993f l or a Ban d function combined with a blackbody an d 
a power law (|Ryde fc Pe’eii [200^ iGuiriec et al.ll201^ . 
The spectrum peaks at photon e nergy Enk which is 
normally comparable to 1 MeV (|Kaneko et al.l 120061 
iColdstein et al.ir2012fl and almost neyer exceeds 10 MeV. 

1.1. Transparent or opaque source? 

A simple interpretation of the obseryed nonthermal 
spectrum would be synchrotron emission from suddenly 
accelerated and quickly cooled electrons in an optically 
thin source. Then the MeV peak of the GRB spec¬ 
trum comes from the peak of the electron energy dis¬ 
tribution ne(7) and the gamma-ray tail comes from the 
high-energy tail in Ueij). This picture has to assume 
that the plasma has no dominant Maxwellian compo¬ 
nent. Fitting the prompt data with synchrotron models 
using a Maxwellian plus power-law electron distribution 
either requires the Maxwellian to be significantly weaker 
than the non-therm al component ([Baring fc Brabyll2004 
iBurgess et al.l 2Ml^ ■ or at most comparable in energy 
( Burgess et al.ll20T€) . The peak of the GRB spectrum 
is then associated with a cutoff of the dominant non- 
thermal electron distribution at low energies. The origin 
of such a special distribution is unknown. Simulations 
of possible heating mechanisms — shocks and magnetic 
reconnection — give different Ueif). Shock heating re¬ 
sults in a dominant Maxwellian component with only 
a s mall fraction of particles po pulating the high-energy 
tail (|Sironi fc SDitkovskvll2011ll . and reconnection gives 
a broad flat distrib ution with no cutoff at low energies 
(|Kagan et al.ll201^ . 

Regardless of the heating mechanism, the synchrotron 
model faces a problem when it is compared with obser¬ 
vations. The synchrotron spectrum has a rather broad 


peak, even when the source is assumed to have a uni¬ 
form magnetic field B and the power-law electron dis¬ 
tribution is assumed to have a sharp cutoff at low en¬ 
ergies. Attempts to reproduce GRB spectra with such 
idealized mo dels gave acceptable fits for a small num- 


ber of bursts (iBarini 

g & Brabvll2004 IBurgess et al.ll201lL 

12014 iPreece et al.1 1 

2014f). In a more realistic model, 


several factors inev itably broaden the synchrotron peak 
(|Beloborodovl [2M^ hereafter B13), which makes it in- 
compatible with observed spectra in virtuall y all GRBs 
(|Axdsson fc Borgonovol[2013 : lYu et al.ll2015f) . 

This problem is illustrated in Figure [I] where a syn¬ 
chrotron spectrum is compared with the Band fit of a 
bright burst with a typical spectrum, GRB 990123. In 
the realistic fast-cooling regime, the minimum width of 
the synchrotron peak at half maximum exceeds 1.5 or¬ 
ders in photon energy, even with a single value of B in 
the emitting source and the narrow injected electron dis¬ 
tribution (Maxwellian). It is significantly broader than 
the Band fit to the time-ave rage of the variable GRB 
spectrum ([Briggs et al.lll99^ . Uncertainties in the ob¬ 
served spectrum due to the detector response and limited 
photon statistics allow some room for stretching the mea¬ 
sured peak width, however can hardly be made consis¬ 
tent with the synchrotron model, especially when the in¬ 
evitable broadening due to variable magnetic field, elec¬ 
tron injection, and Doppler factor is taken into account. 

The sharp MeV spectral peak provides strong evidence 
for thermalization of radiation at early, opaque stages 
of the GRB explosion. The inheritance of the spectral 
peak from an initial thermalization stage is also sup¬ 
ported by the observed distribution of i?pk, which cuts 
off above ~ 10 MeV in agreement with a theoretical max¬ 
imum (B13). In addition, many GRBs show hard spec¬ 
tral slopes violating th e “line of death” fo r optically thin 
synch rotron emission ([Preece et al.l 120001 : iKaneko et al.l 
l2006f) and suggesting an opaque source. 

Accepting that the MeV peak of GRB spectrum forms 
inside an opaque jet leads to so-called “photospheric” 
emission models. The burst radiation is released where 
the jet becomes sufficiently transparent to scattering, and 
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Figure 1. Synchrotron spectra from an optically thin spheri¬ 
cal shell with three electron distributions: mono-energetic (dotted 
line), Maxwellian (short-dashed line), and fast-cooling Maxwellian 
(long-dashed line). For comparison the Band fit of GRB 990123 is 
shown by the red thick curve. 


its spectrum is mainly shaped by subphotospheric ra¬ 
diative processes. Several versions of the p hotospheric 
model have been proposed over the years (iThp imjspd 


19941 lEichler & LevinsonI 

20001 IMeszaros & ReesI 120001 

Giannios & Soruit 12007: 

Beloborodov 20^ hereafter 


B10: lLevinsoDll2012tlThomDSon fc Gillll2014D . All share a 


key feature: the jet is dissipative, i.e. significantly heated 
as it propagates away from the central engine. This 
heating modifies the emitted photospheric radiation from 
simple blackbody emission. The resulting spectrum was 
shown to have a nonthermal shape that clo sely resem¬ 
bles t he phenomenol ogical Band function (iPe’er et al.l 
1200(31 : iGiannio j [20081: BIO: iVurm et al.ll201lL hereafter 
VI1). It was proposed that the dissipative photosphere 
mod el provides a good description to the observed spec¬ 
tra (jBvde et al.l 120111 ) and needs to be carefully tested 
against observations. 


Additional evidence for dissipation at small radii is pro¬ 
vided by the observed photon number emitted in GRBs. 
In many GRBs, the central engine is unable to pro¬ 
vide the observed photons, so additional photons must 
be produced in the expanding jet. Photon production 
is a direct c onsequence of diss ipation at large optical 


depths (B13: IVurm et ahl 

20131 hereafter V13: see also 

lEichler & LevinsorJ 120001 

Thomnson et ahl 120071). Ob- 


servations also require that dissipation continues at least 
to the photospheric radius, so that the released spec¬ 
trum has a nonthermal shape. Therefore, in this paper 
we consider outflows which remain dissipative across a 
broad range of distances from the central engine, start¬ 
ing from the region inside the progenitor and extending 
to the jet photosphere and beyond. 

As long as baryons dominate the plasma inertia, dissi¬ 
pation of internal motions may be expected to heat the 
ions (and neutrons) in the first place. Efficient dissipa¬ 
tion should give a typical energy of ~ 1 GeV per nucleon 
(its rest mass) in a relativistic jet. Baryons themselves 
do not emit significant radiation, because of their large 
mass-to-charge ratio, however their energy can be passed 
to the electrons and radiated in the following ways: 

1. Goulomb collisions gradually pass energy from the 
hot ions to the thermalized electron/positron pop¬ 
ulation (which is kept much colder by efficient ra¬ 
diative cooling). 

2. Inelastic (pion-producing) nuclear collisions gener¬ 
ate a non-thermal population with Lorentz fac¬ 
tor 7 inj ~ vriTrlme ~ 300. 

3. Plasma motion through a radiation-mediated shock 
has a steep velocity gradient, leading to bulk 
Gomptonization of photons by the electrons. Bulk 
Gomptonization could also result from plasma tur¬ 
bulence on small scales. 

4. Gollective plasma processes in a collisionless shock 
suddenly heat the electrons to an ultra-relativistic 
temperature. 


1.2. Internal dissipation 

Four dissipation mechanisms have been proposed 
as sources of GRB e mission: collisionless shocks 
l|Rees fc Meszaro ^ I1994D , da mping of Alfven wave 
turbulence (iThomosonl 19941) . magnetic reconnection 
(|Drenkhahn fc SpruitI I2002D . and neutron collisions 
(BIO). Magnetic field and internal bulk motions provide 
the energy reservoirs available for dissipation. 

The presence of strong internal motions in the jet is 
indicated by the observed variability of GRB radiation. 
The central engine of the explosion is likely unsteady, 
and additional variability is induced as the jet burrows its 
way through the progenitor s tar and the cocoon p roduced 
by the jet-star interaction (iLazzati et al.l 1200^ . This 
leads to multiple internal and recollimation shocks, which 
keep the jet hot and relatively slow when it emerges from 
the stellar progenitor. Thus, shock heating is expected 
to occur in an extended range of radii and in an extended 
range of timescales, which is consisten t with the observed 
broad power spectrum of variability (|Beloborodov et ahl 
120001 : iMorsonv et al.l[20T^ . 


All these proc esses are important in sub-photospheric in¬ 
ternal shocks (jBeloborodovll201^ . Electrons can also be 
directly heated by magnetic reconnection, tapping into 
the magnetic energy carried by the jet. 

1.3. Evolution of radiation in the expanding jet 

The energized electrons rapidly lose their energy to ra¬ 
diation via inverse Compton (IC) scattering, synchrotron 
emission, and (at extremely high optical depths) through 
double Compton scattering and bremsstrahlung. The 
produced photons are redistributed in energy by Comp¬ 
ton scattering and form the spectrum that eventually 
escapes at the Thomson photosphere A* where the scat¬ 
tering optical depth tt drops below unity. 

Three relevant regions in the jet were described in B13: 

1. The Planck zone (r < 10^° cm, tt > 10®): the 
density of the jet is sufficiently high to maintain 
blackbody radiation in detailed equilibrium with 
the thermalized plasma. 
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2. The Wien zone (tt > 10^): the dissipated heat is 
thermalized into a Bose-Einstein photon distribu¬ 
tion with a finite chemical potential. The number 
of photons accumulated in the Wien peak attains 
its final value near the Wien radius, beyond which 
Comptonization is unable to bring new generated 
photons to the spectral peak. 

3. Unsaturated Comptonization zone (tt < 10^): 
heating maintains a Compton parameter y ~ 1. 
The final non-thermal shape of the spectrum is pro¬ 
duced in this region. 

The CRB radiation emerging at the photosphere carries 
information about the entire expansion history of the jet. 
Thus, the observed spectrum may be used to reconstruct 
dissipative processes hidden in the opaque region behind 
the photosphere. 

It is best to view the formation of the observed spec¬ 
trum as a problem of radiative trans fer in an ultra- 
relativistic outflow (|Beloborodovl l20lTl hereafter Bll). 
The problem must be solved consistently with the flow 
dynamics and heating. Boundary conditions should be 
set sufficiently close to the central engine where the hot 
and dense plasma is capable of producing and thermal- 
izing photons, so that these processes can be explicitly 
followed by the transfer simulation. Radiative transfer 
from the inner region to and beyond the photosphere 
determines the radiation spectrum released by the jet. 

1.4. This paper 

Our main goal in this paper is the development of a 
global radiative transfer model for heated jets and the 
application of the model to observed GRBs. The main 
novel feature of our simulations is the explicit inclusion 
of photon production reactions in the transfer problem. 
This requires one to start the simulation at a small ra¬ 
dius, extremely deep below the photosphere. 

This paper focuses on jets whose energy is domi¬ 
nated by matter and radiation, and a modest magne¬ 
tization parameter eb < 1 is assumed. We defer our 
study of magnetically-dominated jets to a future work; 
a recent discussion of radiati ve processes in magnet i¬ 
cally dominated je t s is fo und i niGill fc ThompsonI (1201411 : 
iThompson k, Gilll (|2014ll and iBeeue &: Pe’ei (|2015ll . 

To keep the number of parameters to a minimum we 
employ a simple model of continuous internal dissipation 
approximating the average heating rates (thermal and 
non-thermal) as power laws of radius. This may be a 
crude approximation to jet heating by multiple internal 
shocks or reconnection, however it allows us to study the 
global picture of the evolution of radiative processes with 
radius. As will be demonstrated below, the formation of 
photospheric radiation extends over several decades in 
radius, and the transfer simulation allows us to study all 
stages of this process. The small number of parameters 
in the model makes it useful for fitting the data. 

The simulations presented in this paper are per¬ 
formed with an improved version of the kinetic code 
develo ped by IVurm fc Poutaneiil (120090 and IVurm et’aTI 
(|201lO . It follows in detail all relevant radiative pro¬ 
cesses expected in a heated jet, thermal and nonther- 
mal. The thermal processes include Gomptonization of 
photons by the thermal plasma, induced Gompton scat¬ 
tering at low frequencies, cyclotron emission/absorption. 


bremsstrahlung, and double Gompton scattering. The 
rates of all these processes are well defined and accu¬ 
rately calculated for a given density, temperature, and 
magnetic field. 

The injection of relativistic leptons by nuclear colli¬ 
sions or by collisionless heating (Beloborodov 2016) leads 
to synchrotron emission inside the jet, which is an impor¬ 
tant source of photons. Relativistic leptons also generate 
an IG cascade; it is calculated in detail using exact 
cross sections. The density of pairs is governed by the 
rates of their creation and annihilation; accurate calcula¬ 
tion of density is essential as it can strongly dominate 
over the electron-ion plasma density. 

The paper is organized as follows. Section [2] describes 
the setup of the transfer problem and the physical pro¬ 
cesses involved. Sample models are presented in Section[3] 
where we explore how the emitted spectrum is formed as 
the jet expands over several decades in radius, how the 
emerging radiation is influenced by the jet magnetiza¬ 
tion, acceleration history etc. In Section |31 we apply the 
model to three bright bursts with high-quality spectral 
data: GRB 990123, GRB 090902B and GRB 130427A. 
The results are discussed in Section [SJ 

2. SETUP OF THE TRANSFER SIMULATION 
2.1. Dynamics of the jet 

We assume that the jet collimation occurs inside a 
characteristic radius Rc and its subsequent expansion at 
r > Rc may be approximated as conical. In our models, 
Rc is typically chosen around 10^^ cm; our simulations 
start at this radius. The jet is still accelerating at this 
stage, and its Lorentz factor saturates at a much larger 
radius; the acceleration is self-consistently calculated as 
described below. 

The jet variability is viewed as internal turbulent mo¬ 
tions AT superimposed on a steady flow with Lorentz 
factor r(r). At any radius r, the internal motions can be 
divided into two parts: large-scale “frozen” fluctuations 
(length-scale larger than the local “horizon” — the size 
of the casually connected region) and small-scale evolv¬ 
ing fluctuations. This decomposition is similar to the 
analysis of primordial perturbations in cosmology. The 
small-scale fluctuations are dissipative, and they are as¬ 
sumed in our radiative transfer problem to be the source 
of heat and nonthermal particles. The jet heating must 
be unsteady; however, we simplify the simulation by us¬ 
ing a smoothed, averaged heating rate. When calculating 
radiative transfer, we average out AT and view L as a 
single-valued function of radius r (or the comoving time 
of the steady background flow). Then radiative transfer 
can be modeled in a steady-state approximation, which 
significantly simplifies the calculations (see Bll for dis¬ 
cussion of this approximation). 

The turbulent motions AT are only viewed as a reser¬ 
voir of energy available for dissipation. Then the total 
jet luminosity may be decomposed as 

T — Lp\ Trad “b -^B “b A^urb “ COnst. (1) 

Here Lp\(r) is the steady component of the plasma en¬ 
ergy flow rate separated from the turbulent component 
Aturb(?') (the reservoir); Arad('c) and Lsir) are the energy 
flows carried by radiation and magnetic field. As the jet 
propagates, energy can be redistributed between the dif- 
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ferent components in Equation (III)- The dissipated part 
of the turbulence reservoir gets converted to the plasma 
internal energy, which is immediately transferred to radi¬ 
ation via rapid cooling; thus effectively there is a gradual 
transformation Lturb Trad- In the simulations pre¬ 
sented below, we keep Lb = constant for simplicity; in 
reality, the magnetic component may have a ’reducible’ 
component available for dissipation, e.g. via reconnec¬ 
tion. However, when eb = Lb/L 1 the dissipation of 
magnetic fields cannot serve as the main source of GRB 
radiation. 

Radiation can efficiently accelerate the jet at early 
stages, when it dominates the jet energy. Thus we also 
expect the gradual transformation Trad Tpi. The dy¬ 
namical equation for the jet Lorentz factor E is given by 
(see Bll), 


dr „ Airh 

— = (TT^i -^ 

dr rripC^ 


( 2 ) 


where 47r/i is the radiation flux measured in the rest- 
frame of the jet (Ji is the first moment of the local radi¬ 
ation intensity), and Z± is the number of electrons and 
positrons per proton. Equation ([2]) assumes Tpi Ri TMc^; 
this is a good approximation as long as the ion temper¬ 
ature is non-relativistic. 

We use Equation (I2|) in our simulations to calculate the 
self-consistent r(r). At very large optical depths, where 
radiation is nearly isotropic and Ii is small, an alterna¬ 
tive form of the dynamical equation can be use d, wh ich 
does not depend explicitly on /i, see Equation (IA13I1 in 
Appendix [A) 


2.2. Kinetic equations 

At optical depths tt > 100 the radiation field is nearly 
isotropic in the jet rest frame, and the bulk of leptons are 
kept at a non-relativistic temperature. The evolution of 
radiation under such conditions is well described by the 
Kompaneets equation. For a conical accelerating outflow 
it takes the form 




x'^ dx 


^ I „ dh _ 3 — g 

ttx { 0e-^+n + n^ j - — 


X^fi 


(3) 


where n is the photon occupation number, x = hiy/rUeC^ 
is the dimensionless photon energy and 0g = kBTg/mgC^ 
is the dimensionless electron temperature. The last term 
in square brackets accounts for adiabatic cooling of radi¬ 
ation, where 


to an overestimation of the photon number in the Wien 
peak. 

The deviation from isotropy develops at optical depths 
Tt ~ a few tens (Bll), roughly near the Wien radius. 
At larger radii one has to s olve the full an gle-dependent 
radiative transfer equation (jMihalaslIlQSDl : B11)I3 


1 ^ V 

f) T f) 

+ (1 + g)(l - 9L)-gJ^ “ ^ 


(5) 


where I^, is the comoving specific intensity, g = cos 0, 0 
is the (comoving frame) angle relative to the radial direc¬ 
tion, and ji, is the emissivity. This equation accurately 
describes the relativistic radiative transfer in three di¬ 
mensions for spherically symmetric outflows with T > 10. 

The evolution of the electron/positron distribution is 
described by the kinetic equation 


1 (9 v 

^ " CK±n±(p)) 


d ( r 
dp \ cT 


1 d 

pn±(p) - - — [Dn±{p)] 


3-g 


pn±{p) 


( 6 ) 


where p = \ is the electron/positron momentum 

in units of rUeC, and n±{p) is the electron/positron den¬ 
sity per unit momentum interval. The terms p and D 
account for heating/cooling and diffusion in momentum 
space due to radiative processes and Coulomb collisions. 
The emission term j± includes the injection of new pairs 
due to non-thermal dissipation and photon-photon ab¬ 
sorption 7 + 7 —>■ e^ + e", and CK±n±{p) is the pair anni¬ 
hilation rate. All these terms are accura tely calculated as 
descri bed in IVurm fc Poutaneiil (1200911 and I Vurm et'ffi] 
(1201111 . The last term in curly brackets describes adia¬ 
batic cooling. 

In the simulation the photon equation is solved to¬ 
gether with the kinetic equations for electrons/positrons 
and the dynamical equation ([2]) for r(r). The photon and 
lepton equations are coupled via the interaction terms jn 
{jv), Ku, j±, K±, p and D. The switch from the Kom¬ 
paneets equation ([3]) to the radiative transfer equation 
(O is made at tt = 100; at this stage either of these 
equations is sufficiently accurate, allowing for a smooth 
transition. The evolution is followed to tt «C 1. The 
simulation is typically stopped at tt ^ 0.1 where the 
obtained radiation field is transformed to the observer 
frame and assumed to escape. 


All quantities except r and T are measured in the instan¬ 
taneous comoving frame of the jet. The opacity k^, and 
the source term j^, account for all relevant radiative pro¬ 
cesses (see Section [23|l except thermal Comptonization, 
which is described by the Kompaneets equation itself. 
Note that it is essential to include the induced term 
in the Kompaneets equation to accurately account for the 
number of synchrotron photons that can be upscattered 
to the Wien peak (V13); neglecting this term would lead 


2.3. Radiative processes and photon generation 

Radiative processes in the expanding jet determine the 
observed GRB spectrum and require accurate calcula¬ 
tions at all radii, including very opaque regions far be¬ 
low the photosphere. In particular, the position of the 

^ In contrast to the Kompaneets equation, Equation © neglects 
induced downscattering, which affects the rate of photon upscatter- 
ing from low energies toward the Wien peak. This loss of accuracy 
is acceptable, because our simulation switches to Equation © out¬ 
side the Wien zone where upscattering has stopped populating the 
Wien peak with new photons. 
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spectral peak is controlled by photon generation below 
the Wien radius, where photons are produced at low en¬ 
ergies and up-scattered to the peak via saturated Comp- 
tonization. The photon generation mechanisms can be 
categorized a follows (B13, V13): 

1. Thermal; photons are generated by Maxwellian 
leptons. Thermal radiative processes operate ef¬ 
ficiently only in the Planck zone (tt' > 10®). 
These processes include double Compton scat¬ 
tering, bremsstrahlung (less efficient than double 
Compton unless the jet is strongly clumped), and 
cyclotron emission. The efficiency of cyclotron 
emission can be similar to that of double Comp¬ 
ton scattering if the jet is strongly magnetized. 

2. Nonthermal: photons are generated by relativistic 
nonthermal electrons/positrons. The main radia¬ 
tive process is synchrotron emission. It continues 
to operate outside the Planck zone, and the num¬ 
ber of produced photons typically increases with 
radius, as long as the non-thermal dissipation chan¬ 
nel remains active. 

The rates of the thermal processes are given in B13 and 
V13. The implementation of nonthermal synchrotron 
emission is similar to that in VI3, with one significant 
extension: we are considering jets with moderate mag¬ 
netization £b ^ 10“® — 10“® where pair-photon cascades 
produce a larger number of secondary nonthermal par¬ 
ticles (see below). In contrast, VI3 considered magne¬ 
tizations close to equipartition where pair cascades are 
suppressed. 


2.4. Heating rate 

The total (thermal and non-thermal) heating rate sat¬ 
isfies 

dTji _ dZ/tB ^ dZ/nth _ dZ/turb 

d In r din?’ d In r d In r 

The simulations presented below assume that the frac¬ 
tional luminosities dissipated via the thermal and non¬ 
thermal channels follow a power law, 


^nth — j- 


I 

( r 

(8) 

L dinr 



/ r \ 

(9) 

L dinr 

= e0.nth(^^j • 


The corresponding heating and injection rates in Equa¬ 
tion m are 


and 


7 . _ 7_ SthL _ 

p p 47r?7ieC^r®r(r?_-|-7r_|_) ’ 


( 10 ) 


_ ^nthL 5{j) Pinj) 

47r?7leC^r®r 27 inj 


( 11 ) 


where n- and n+ are the electron and positron densities, 
respectively, pinj is the dimensionless momentum of the 
injected pairs, and 7inj = (I + i® their Lorentz 

factor. 


2.5. Parameters of the model 

The GRB models presented below have the following 
set of parameters: 

(1) Total jet luminosity L (isotropic equivalent). 

(2) Energy per baryon ry = L/Mc^. 

(3) Magnetization eb = LbIL. 

(4) The jet Lorentz factor at the radius Rc where the 
simulation starts; Rc is chosen near the radius of the jet 
breakout from the progenitor star. 

(5) The jet energy per photon, i.e. the initial Epk. In our 
models the jet is initially photon starved, i.e. we choose 
a high Epk ~ 10 MeV. This is done to demonstrate that 
even in weakly magnetized jets synchrotron emission can 
efhciently produce photons. Then the predicted GRB 
spectrum weakly depends on the choice of initial Zipk. 

(6) The efhciencies of thermal and nonthermal heating, 
Eth and Enth- In niost of our models heating is assumed 
to have a flat distribution over Inr (Zth = ^nth = 0 in 
Equations ([8|) and (H))); we will also consider a few vari¬ 
ations in the heating history of the jet if this is required 
to fit the observed GRB spectra. 

The non-thermal lepton injection energy is fixed at 
7inj = 300 ~ ?7i,r/nT.e, as expected for collisional heat¬ 
ing; its exact value weakly affects the results. 

3. SAMPLE MODELS 
3.1. Overview of the spectral evolution 

Figure [2] illustrates the evolution of the photon spec¬ 
trum within a dissipative jet as the radiation is carried 
from highly opaque regions to transparency. There are 
two main stages of the spectral evolution: 

(1) Generation of photons below the Wien radius (tt > 
10^) and their Gomptonization to the Wien peak. The 
growth of photon number in the Wien peak results in 
its shift to lower energies. The high-energy tail is sup¬ 
pressed at r < i?w for two reasons. First, the electrons 
are kept approximately at the Gompton temperature, i.e. 

ft! Epk/r. Thus thermal Gomptonization is hardly 
capable of populating the photon spectrum above Epk, 
leading to an exponential cutoff at Epk. Secondly, the 
nonthermal radiation from IG scattering by relativistic 
pairs does not survive at tt » 1 — it is reprocessed 
to lower energies via the pair-photon cascade and down¬ 
scattering by the thermal pairs. 

(2) Broadening of the spectrum by unsaturated Gomp¬ 
tonization outside the Wien radius, leading to a nonther¬ 
mal shape of the spectrum. The broadening first affects 
the low-energy slope a, which begins to soften near the 
Wien radius and attains its final value near the photo¬ 
sphere. The resulting photon index a ^ — 1 depends on 
the heating history at r > R-^, in particular on the non¬ 
thermal dissipation channel, as the unsaturated Gomp¬ 
tonization of the fresh low-energy synchrotron photons 
governs the softening of a. Without nonthermal dissipa¬ 
tion the low-energy spectral slope would be much harder. 

The photospheric spectrum also broadens at high en¬ 
ergies. The high-energy tail is mostly built around and 
above the photosphere, due to two effects, thermal and 
nonthermal. Outside the Wien radius the electron tem¬ 
perature rises significantly above Zipk/E (see below) and 
thermal Gomptonization begins to produce photons with 
energies E > Zipk. As the optical depth drops, the 
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E [MeV] 


Figure 2. Evolution of the radiation spectrum carried by the 
jet, from high to low optical depths tt- This sample model has 
the following parameters: L = 2 X 10®^ erg/s, rj = 190, Rc = 
10^^ cm, 70 = 300, eth + Enth = 0.05, eth/£nth = 1, r(i?c) = 20, 
magnetization eb = 10“^. The spectra are measured in the rest 
frame of the central engine and not corrected for a cosmological 
redshift. 

nontliermal high-energy component becomes increasingly 
prominent, especially in weakly magnetized jets. Then 
the overlapping thermal and nonthermal Comptonization 
components together form an extended high-energy spec¬ 
trum which can easily be mistaken for a single emission 
component (see also BIO; Vll). In contrast, in strongly 
magnetized jets with ep 1 , the pair cascade and non¬ 
thermal Comptonization are reduced, leading to a dif¬ 
ferent spectral shape at high energies (see Figure 5 in 
Vll). 

It is worth noting that the overall spectrum has a dis¬ 
tinctly nonthermal shape already well below the photo¬ 
sphere (at Tt ^ 10). In particular, the low-energy slope is 
significantly softer than the Planck (or Wien) spectrum. 
Even if dissipation stopped completely at tt ~ 10, a 
thermal-looking prompt emission would not be expected; 
instead, the emerging spectrum would resemble a cutoff 
power law. 

A thermal-looking GRB spectrum would only be pro¬ 
duced in a special case when all dissipation is conhned to 
very high optical depths. Even in this case, however, the 
transfer effects would soften a from the th ermal slope to 
a ~ 0.4 (BIO; see also lPe’er fc Rvdjl2011[l . 

3.2. cascade and pair loading 

When the jet magnetization is weak, eb < 0.1, the rel¬ 
ativistic pairs injected by the nonthermal dissipation 
channel lose their energy primarily through IC scatter¬ 
ing of photons from the quasi-thermal (Wien) peak. The 
upscattered photons with energies above meC^ in the jet 
frame immediately convert to secondary pairs through 
photon-photon collisions. The created particles also up- 
scatter photons etc., leading to a pair-photon cascade 
that populates the jet with secondary pairs which 
can dominate the plasma. 

The distribution of non-thermal leptons resulting from 
the cascade can be expressed as (cf. Equations (18) and 


(19) in Vll) 


neil) 


4cCTt (Eb + £rad) r 


— (a±+2) 

I 5 


(12) 


where Eb = L ^dLn/dl nr. Eras = L ^d Araa/dlnr, and 
Y ^ 0.1 is the pair yield (ISvenssonlll987f) . The saturated 
cascade turns off at 7 ±,min ~ (weC^P/Epk)^/^ ^ 10 ; elec¬ 
trons at 7 < 7±,min are unable to upscatter photons from 
the Wien peak to MeV energies for further pair genera¬ 
tion. The pair yield of the cascade generated by a pri¬ 
mary particle with Lorentz factor 7 inj may be expressed 
as V = A4s/7inj, where A4s is the pair multiplicity, i.e. 
the number of secondary per primary particle. The 
power-law index a± In Als/(ln 7 inj — ln 7 ±,min) charac¬ 
terizes the additional steepening effect of the cascade 
on the standard cooling distribution ne{'y) oc 7 “^. 

Below the photosphere the pairs are in creation- 
annihilation equilibrium. Their density can be expressed 
as 

Tie — 

where the proton density is 


M 

dirrupcr^r ’ 


(14) 


and the pair loading factor Z+ is give n by (Appendix [BJ 
see also BIO: iThomnson fc Gil]ll20fll 




( leVEnthl^ Wp 
ySrpTMc^ TOe 


1/2 


oc 


(rr)i/2. 


(15) 


The last proportionality is valid if Z± ^ 1, i.e. if pairs 
dominate over the electron-ion plasma. This condition 
is met when the Thomson optical depth associated with 
Up is sufficiently low. 


o'Tnpr ^ 16 TWp YenthL 

r 3 me TMc^ ■ 


(16) 


This condition may not be satisfied in a region deep be¬ 
low the photosphere. The expanding jet becomes increas¬ 
ingly pair dominated in the subphotospheric region, pro¬ 
vided that the nonthermal dissipation channel is active 
(Figure [21). 

Pair annihilation freezes out when the annihilation 
time, Gnn = 16/(3c(Tt?t-±): exceeds the expansion time 
tdyn = f/icT), which occurs at t± = 16/3. The pair¬ 
loading factor at this time is 


£nthVT 

me FMc^ 


(17) 


where F, Enth and Y are taken at the freezeout radius. 
Typically, we find Z± ^ 10 for moderate nonthermal in¬ 
jection rate Enth ~ 0.1. As shown in Appendix [Bland seen 
in Figure m pair creation beyond the freezeout radius can 
increase Z± by a modest factor (at most logarithmically 
in r). 

The pair loading increases the photospheric radius by 
the factor of Z± compared to the pair-free jet. Pair load¬ 
ing is significantly reduced if the jet magnetization eb is 
above a few percent, since synchrotron cooling becomes 
competitive with Compton cooling, inhibiting the pair- 
photon cascade and lowering Y (Vll). 
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Figure 3. Evolution of pair-loading with radius in the simulation 
shown in Figure[2] The vertical dashed line indicates the freeze-out 
radius where t± = 16/3. 


3.3. Building up the photon number in the Wien peak 
via synchrotron emission 

A fraction of the non-thermal lepton energy is con¬ 
verted to low-energy synchrotron radiation. If the mag¬ 
netization is weak, this fraction is small, as the syn¬ 
chrotron losses of high-energy particles are small com¬ 
pared to their IC losses. While the energy budget of syn¬ 
chrotron radiation is small, the number of generated syn¬ 
chrotron photons is substantial, because of their low en¬ 
ergies. Furthermore, these soft photons can gain energy 
from the thermal plasma through Comptonization. Be¬ 
low the Wien radius i?w many of them are Comptonized 
to the Wien peak and eventually dominate the peak, sig¬ 
nificantly increasing its photon number and shifting Apk 
to lower energies. This shift is required by energy con¬ 
servation (energy is shared between more photons) and 
accomplished through Compton cooling of the thermal 
plasma by the synchrotron photons. 

The number of synchrotron photons that can be up- 
scattered to the Wien peak depends on the competition 
of Compton upscattering against synchrotron reabsorp¬ 
tion and induced Compton downscattering (V13). This 
competition defines a critical photon energy Eq (and a 
corresponding Lorentz factor 70 ), above which upscat¬ 
tering dominates over other processes. The calculation of 
7 o, Eq, and the number of synchrotron photons emitted 
above Eq is given in Appendix [Cj 

At r < all synchrotron photons emitted at E > Eq 
end up in the Wien peak. Most of them are emitted near 
Eq by leptons with Lorentz factors 7 ~ 70 (the emissivity 
scales as where p > 2). Therefore, the rate at 

which photons are accumulated in the peak [cm“^ s“^] 
is roughly proportional to the number of leptons with 
7 ^ 70, 


A4 


synch ‘ 


aTUleC^B 

3eh 


7o’^e(7o)- 


(18) 


The number of photons accumulated in the Wien peak is 
hnalized near The resulting flow of photon number 
through the sphere AttR^ (isotropic equivalent) is given 


by 


/ Rw 

Msynch{r)d\nr. (19) 

The photon number carried in the spectral peak is weakly 
changed outside the Wien radius; this number is eventu¬ 
ally released at the photosphere. 

The secondary pairs affect photon production by 
increasing both and 70 ; the net effect is positive, i.e. 
the pair cascade enhances photon generation. Pair load¬ 
ing also increases i?w by increasing the flow opacity; this 
allows a longer time for the accumulation of photons in 
the Wien peak. 

In contrast to thermal bremsstrahlung or double 
Compton effect, the synchrotron photon production in¬ 
creases with radius and the integral in Equation (I19|) 
peaks near i?w This behavior was seen for pair-free 
outflows in V13 and remains true in the presence of pair 
cascades (see Appendix [Cl for analytical estimates). Our 
numerical simulations confirm that a substantial fraction 
of the photons accumulated in the spectral peak origi¬ 
nates near the outer boundary of the Wien zone. This 
can be seen in Figured the decrease of Epk due to con¬ 
tinuing photon supply to the peak ends around tt < 100 , 
which approximately corresponds to i?w- 
As long as the nonthermal dissipation channel remains 
active, the production of synchrotron photons continues 
beyond the Wien radius. This results in the soft “ex¬ 
cess” seen below a fewxlO keV in the emitted spectrum 
(Figure m see also VI1). More importantly, unsaturated 
Comptonization of these photons plays a key role in de¬ 
termining the spectral slope below the peak, which is 
discussed below. 


3.4. Comptonization 

The saturated Comptonization at r < i?w maintains a 
quasi-equilibrium between radiation and thermal plasma. 
The photon spectrum around the peak has a Wien shape, 
with small excesses at lower and higher energies due to 
synchrotron and non-thermal IC emission, respectively. 
At r > i?w, the thermal Comptonization gradually 
switches to the unsaturated regime, where the Comp¬ 
ton y-parameter remains close to unity, as long as the 
heating operates. 

The y-parameter may be evaluated by considering the 
evolution of the radiation luminosity at tt ^ 1 (see Ap¬ 
pendix!^, 


dTrad 2y 4 dLk 

d In r 3 3 d In r 


( 20 ) 


where the first and second terms on the right hand side 
account for adiabatic cooling and dissipation, respec¬ 
tively. Assuming g = 1 — d In F/d In r Ri const and the 
scaling 


dLk 




dlnr 

one can solve Equation (1^ for Trad, 

4 dLh 


Trad — 


3fch -b 2y dlnr 


( 21 ) 


( 22 ) 
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The (thermal) heating of electrons is balanced by Comp¬ 
ton losses, 

^ = (23) 

where y = drrSe and yc = 4tt0c are the Compton 
parameters corresponding to the electron and Comp¬ 
ton temperatures, respectively. Substituting the solution 
into the heating-cooling balance (1251) . one finds that 
the Compton parameter relaxes to 

y-yc = kh + Y- (24) 

Thus y—yc remains close to unity if the heating rate does 
not decline significantly faster than dL\i/d\nr oc r°. At 
the Wien radius y ~ yc ~ a few and 0e ~ ^c- At larger 
radii the electron temperature 0^ increases above the 
Compton temperature 9c-, as Comptonization gradually 
switches to the unsaturated regime with y Ri /ch + 2y/3. 

The region with y ^ 1 can extend over several expan¬ 
sion times. Therefore, the effective net Compton param¬ 
eter can be substantially larger than unity and result in 
significant redistribution of photons, in particular above 
the spectral peak. The nonthermal high-energy tail is 
built in this regime, with an energy content comparable 
to that of the radiation near the peak. In addition to 
thermal Comptonization, the pair photon cascade initi¬ 
ated by the injected non-thermal leptons can extend the 
non-thermal tail to the GeV range. 

The final low-energy slope of the photospheric spec¬ 
trum is mainly formed at radii r ^ i?w, in the interme¬ 
diate regime with y ~ a few. It is shaped by the combi¬ 
nation of photon Comptonization from low energies to¬ 
ward the peak, diffusion in energy space, and adiabatic 
cooling. 

Some insight into the development of the low-energy 
slope is provided by the quasi-steady solution of the 
Comptonization problem at a given radius in the spec¬ 
tral range Eq < hv < k^Tg. The solution is found by 
setting the second term on the right side of Equation ([3]) 
to zero. Neglecting induced scattering, one findt0 

oc C-b (25) 

where C = const and a = 3 — 4(3 — y)/(3y). In the 
strongly saturated regime y ^ 1, the last term in Equa¬ 
tion (1331) gives the Wien spectrum with a ~ 3. The other 
term oc results from the steady photon flux in the 
energy space toward the Wien peak. 

The relative amplitudes of the two components in 
depend on the generation rate of low-energy photons and 
the number of photons already accumulated in the Wien 
peak. As y decreases, the Wien peak becomes less pro¬ 
nounced relative to the oc component. The spectral 
hardening toward the peak becomes weaker, because a 
decreases from its saturated value. In addition, the peak 
itself is broadened by the competition between Comp¬ 
tonization and adiabatic cooling. 

Our detailed numerical simulations confirm that the 
average spectrum approaches oc below the peak 

^ The actual spectrum at a given r has time to relax to t he s teady 
state solution only if y ^ 1. When J/ ~ a few, Equation II25II only 
shows the qualitative behavior of the spectrum, but does not give 
its exact shape. 



Figure 4. Effect of different Lorentz factors at Rc on the emitted 
spectrum. Other parameters are the same as in Figure [2] 


for a broad range of parameters. The most important 
requirement for this behavior is the existence of a low- 
energy photon source between i?vv and i?*. The spec¬ 
trum /jy oc corresponds to the photon index a = —1, 
which coi ncides with the aver age photon index observed 
in GRBs (jKaneko et al .11200^ . 

3.5. The role of the early acceleration stage and 
magnetization 

Eigure |4] shows the response of the final (observed) 
spectrum to variations in the Lorentz factor at the col- 
limation radius r(i?c). Lowering L at r < i?w implies a 
higher density and increases the efficiency of photon pro¬ 
duction. In addition, is increased, which increases 
the number of synchrotron photons reaching the spectral 
peak (see Equation IC19I1 . Both effects lead to a lower 
T^pk- 

The spectra become more narrow in jets that are still 
significantly accelerating between and i?*, i.e. those 
with a larger disparity between r(i?c) and the final T. 
During the accelerating stage radiation carries a large 
fraction of the jet energy, which makes it less susceptible 
to spectral redistribution/broadening by dissipation. In 
addition, the ratio is smaller if the jet is still 

accelerating outside i?w, leaving less time for broadening 
the spectrum. 

Figure [5] shows the spectra for different jet magne¬ 
tizations. With increasing eb the peak position shifts 
to lower energies, as more synchrotron photons are pro¬ 
duced. This effect is partially offset by the suppression 
of pair loading (see Section [3T31 and VI1), which reduces 
the Wien radius and the number of synchrotron-emitting 
particles. 

The spectral shape exhibits a characteristic behavior 
as the magnetization is increased: more low-energy syn¬ 
chrotron photons tend to make the spectrum softer be¬ 
low the peak; a stronger low-energy “excess” is also pro¬ 
duced. Above the peak the slope becomes steeper as syn¬ 
chrotron losses reduce IC emission from high-energy par¬ 
ticles. Synchrotron losses also inhibit pair cascade, which 
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Figure 5. Effect of different magnetizations. Other parameters 
are the same as in Figure [2] 

weakens and hardens the high-energy spectrum. As a re¬ 
sult, the nonthermal IC emission creates a stronger up¬ 
ward curvature in the spectrum around 10 MeV. 

4. FITS TO PROMPT GRB DATA 

To test our transfer simulations against observations, 
we chose three bright GRBs with good spectral data and 
distinct spectral shapes: GRB 990123, GRB 090902B, 
and GRB 130427A. Good spectral fits were previously 
obtained with phenomenolo gical models combining Band 
function (|Band et al.ll2009ll . a power law, and sometimes 
a Planck component. Our model will pass the test if it 
is able to reproduce the fits. Formal fitting utilizing an 
appropriate goodness of fit statistic as well as instrument 
response matrices is deferred to a future work. 

The results are shown in Figures [Bl[51 The previous 
phenomenological fits are shown by red dashed curves 
(see Table [2] for fit parameters) and our transfer model 
is shown by solid blue curves. The achieved agreement 
demonstrates that the model is consistent with the data 
and provides estimates for the jet parameters for each 
GRB. The parameters are given in Tabled] 

The transfer model also provides a physical interpre¬ 
tation for the previously suggested phenomenological 
components. In particular, the Band component be¬ 
low ~ 10 MeV in all three bursts results from ther¬ 
mal Comptonization of low-energy photons by the heated 
plasma below the photosphere. The high-energy compo¬ 
nent from nonthermal Gomptonization overlaps with the 
thermal Gomptonization component and smoothly ex¬ 
tends it beyond 10 MeV. The associated spectral hard¬ 
ening is more gradual if the thermal heating persists and 
dominates the dissipation above the photosphere, as in 
e.g. GRB 990123. 

The case of GRB 090902B is particularly interesting, 
as it shows strong deviations from the Band function 
at both low and high energies. These deviations were 
previously modeled as an additional power law compo¬ 
nent that extends from the keV band up to GeV ener¬ 
gies. Our transfer simulations show that the soft and 
hard excesses are produced by different radiative pro¬ 
cesses: synchrotron at low energies and IG at high ener- 



Figure 6. Simulated spectrum of GRB 990123 (solid line), and 
a Band fit (dashed line. lBriggs et aLlll999l l. Jet parameters: Mag¬ 
netization £b = 0.018, simulations starting (collimation) radius 
Rc = 3x cm, initial Lorentz factor r(Rc) = 35, Initial number 
of photons per baryon npj,(Rc)/^p(Rc) = 5x 10^, terminal Lorentz 
factor Tf = 500. The heating parameters are eo^th -bso.nth = 0.26, 
£o thAo nth = 4.7, fcth = fcnth = —0.19. The heating proceeds until 

TT = 0.01. 



Figure 7. Simulated spectrum of GRB 090902B (s olid line), 
and a Band+power-law fit (dashed line, lAbdo et all I2009I . bin 
b). Jet parameters: £3 = 0.012, Rc = 3 X 10^*^ cm, r(Rc) = 
70, npii(iic)/np(Rc) = 10®, Tf = 1220; Heating parameters 

^0,th J- S0,nth — 0.12, £o,th/^0,nth “ 0.85, = ^nth — 0.25. 

The heating proceeds until rx = 0.03. 

gies. However, both are emitted by the same nonthermal 
population, a result of strong nonthermal heating of 
the jet around the photosphere. The simplest mechanism 
of nonthermal heating is the decay of pions produced by 
inelastic nuclear collisions (BIO), which must operate in 
GRB jets unless they are strongly dominated by mag¬ 
netic fields. We find that GRB 090902B has ep ~ 0.012. 
This magnetization explains both the soft and hard ex¬ 
cess. Remarkably, the same value of Eb is required by 
the shape of the Band component observed in this burst. 
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Table 1 

Jet parameters for the fitted GRBs. 


GRB 

L" 

[10“ erg s-'] 

r(Rc)'' 

Tf" 



£0,h* 

^0,th/^0,nth^ 


^nth 

R.' 

[lo'" cm] 

990123 

0.44 

35 

500 

730 

0.018 

0.12 

4.7 

-0.19" 

-0.19" 

1 

090902B 

2.3 

70 

1220 

1740 

0.012 

0.12 

0.85 

-0.25' 

-0.25' 

1 

130427A 

0.43 

64 

400 

506 

0.050 

0.18 

1.4 

0.04”" 

0.04 -2”" 

2 


^Total jet luminosity (isotropic equivalent). 

'^Jet Lorentz factor at Rc. 

‘^Final Lorentz factor Ff achieved by the jet at large radii (calculated self-consistently from the model). 
"^Jet energy per unit rest mass rj. 

® Magnetization ^b- 

fTotal heating rate eo,h = eo,th + eo.nth- 
^Ratio of thermal to nonthermal heating eg th/^o nth- 
''Radial scaling index of thermal dissipation fcth- 
'Radial scaling index of nonthermal dissipation fcnth- 

jphotospheric radius H* (not a free parameter, calculated self-consistently from the model). 
'^Dissipation proceeds until tt = 0.01. 

'Dissipation proceeds until rx = 0.03. 

■"The nonthermal dissipation law changes at Tp = 30. Dissipation proceeds until rx = 2. 


Table 2 

Parameters of the phenomenological fits from the literature. 


GRB 

-Fpk,obs 

[keV] 

at 


T 

1 pi 

T-To" 

[s] 

Instrument 

Reference 

990123 

720 

-0.6 

-3.11 

- 

12 - 45 

CGRO/BATSE 

Bribes et al. ('1999') 

090902B 

908 

0.07 

-3.9 

-1.94 

4.6 - 9.6 

Fermi 

Abdo et al. (2009). bin b 

130427A 

1028 

-0.958 

-4.17 

- 

0-19 

Konus/ Wind 

Golenetskii et al. (2013) 


"Spectral peak position Spk.obs in ths observer frame. 
''Low-energy slope a. 

"High-energy slope /3. 

'ipower-law index in Band-|-PL fit. 

"Time interval of the fitted data relative to the trigger. 



Figure 8. Simulated s pectrum of GRB 130427 A (solid line), and 
a Band fit (dashed line. IGolenetskii et al.ir201,3ll . Jet parameters: 
SB = 0.05, i?c = 3 X lO" cm, r(i?c) = 64, nph(Rc)/n.p(Rc) = 
1.2 X lO"*, Pf = 400; Heating parameters eo.th + ^o.nth = 0.18, 
£o,th/£o nth = l'4i ^th = ^nth = 0.04; the non-thermal dissipation 
law changes into fcnth = ~2 at Tp = 30, the heating proceeds until 
Tx = 2. The black triangle corresponds to the average flux above 
100 MeV observed by Fermi/LAT during the m ost intense phase 
of the prompt emission HAckermann et ffll2014l) . 

The distinct spectral shapes in the three cases are 


mostly the result of different heating histories in the jet, 
as well as different partitioning of the dissipated energy 
between thermal and non-thermal channels. The ab¬ 
sence of an extra component in GRB 990123 and GRB 
130427A implies that non-thermal heating is weak in the 
optically thin parts of the jet. Furthermore, the steep, 
almost cutoff-like appearance of the spectrum above the 
peak in GRB 130427A suggests that most of the dissi¬ 
pation (thermal and nonthermal) is confined to regions 
below the photosphere, as heating at larger radii would 
tend to flatten the high-energy spectrum. This last point 
is quite general and can be applied to other bursts whose 
spectra have a cutoff power-law shape. In contrast, the 
relatively soft low-energy slope of GRB 130427A a Ri —1 
indicates that both dissipation channels have to be effi¬ 
cient at optical depths tt 10 below the photosphere, 
where the spectral slope a is formed. The softening of 
a is also helped by the strong magnetization £b ~ 0.05, 
which increases the supply of soft synchrotron photons. 

Note that all three spectra peak just above 1 MeV, 
as is typical for bright GRBs. The transfer model nat¬ 
urally explains the peak position. It is regulated by the 
photon number in the Wien zone, which is calculated 
self-consistently with no fine tuning. 

5. DISCUSSION 

In this paper we studied the production of radiation 
by opaque, dissipative, relativistic jets. We used radia¬ 
tive transfer simulations to study the main features of 
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the emission mechanism and to compare the theoretical 
spectra with observations. We presented a detailed phys¬ 
ical model that fits the observed GRB spectra and thus 
allows us to measure the main parameters of the jet. 

We used a somewhat simplified description of energy 
dissipation in the jet, with continual electron heating of 
two types, thermal and nonthermal. It is mainly moti¬ 
vated by the minimal heating model of BIO where elec¬ 
tron/positron plasma receives energy from hot baryons 
through collisional processes, which are well defined and 
can be calculated from first principles. However, our sim¬ 
ple parametrization of electron heating (power law scal¬ 
ing with radius, with a given slope and normalization) 
may also accommodate other heating models. 

Another simplification adopted in this paper is the 
steady-state approximation to radiative transfer (see Bll 
for a discussion of this approximation). The theoretical 
spectra calculated in this work represent the average over 
the causal timescale near the photosphere (ft* ^ i?*/r^c, 
which is typically very short in the presented models, 
comparable to a few ms. The spectrum can evolve on 
timescales St > St^,, as observed in GRBs. We believe 
that this evolution is controlled by the noisy process of 
jet formation and mass loading near the central engine, 
which is difficult to predict. 

The main advantage of our transfer simulations is that 
they carefully include relevant radiative processes, so 
that the radiation spectrum is accurately calculated. An¬ 
other advantage is that our dissipative jet model has a 
moderate number of parameters (seven for an unbroken 
dissipation profile and nine if a break is required by the 
data). Therefore, the transfer model can be used to fit 
observations. Unlike fits by phenomenological functions, 
fits by the physical model directly provide estimates for 
the main physical parameters of the GRB jets. 

We find that a continuously heated and moderately 
magnetized jet naturally produces a Band-type spec¬ 
trum, with spectral slopes a, /3, and peak position Apt 
consistent with observations. The bulk of the observed 
photons in a typical GRB do not originate from the cen¬ 
tral engine or the vicinity of the Thomson photosphere; 
instead they are produced in the opaque jet at optical 
depths tt > 100. The broadening of the spectrum into 
the final non-thermal shape takes place between tt ~ 10 
and ^ 0.1. 

A typical burst with L Ri 10®^ erg/s, a canonical Band 
spectrum, and Epk ^ 1 MeV is reproduced by the model 
if: 

1. The jet magnetization is in the range 10“^ < 
£b ^ 0.1. Very weak magnetization increases Apk 
by suppressing synchrotron emission; strong (near 
equipartition) magnetization softens the spectrum 
both below and above the peak, and generates a 
prominent soft “excess” below a few tens of keV. 

2. The jet Lorentz factor T{Rc) < 100 at radii com¬ 
parable to that of the stellar progenitor. Low Apk 
bursts are more “photon-rich” and require consid¬ 
erably lower r(i?c) for more efficient photon pro¬ 
duction. 

3. Heating has a nonthermal component that injects 
relativistic leptons into the jet. The absence of non¬ 


thermal leptons would lead to hard spectra with 
high Upk, due to the lack of synchrotron emission. 

Our simulations demonstrate that nonthermal heating 
in the sub-photospheric region loads the jet with pairs 
via pair-photon cascades. In weakly magnetized jets the 
pairs outnumber protons by a factor comparable to 10, 
leading to the increase of the average photospheric radius 
by a similar factor. The jets remain forever dominated by 
pairs, which will affect afterglow emission produced 
by the reverse shock when the jet is decelerated by the ex¬ 
ternal medium. The reverse shock emission might serve 
as a probe of pair loading. 

5.1. Spectral peak position 

Transfer simulations allow one to study the physical 
conditions in the jet, how the prompt radiation spectrum 
is formed, what controls the observed spectral index, the 
position of the peak etc. 

The evolution of the radiation spectrum in the expand¬ 
ing jet takes place in two stages delineated by the Wien 
radius i?w where Comptonization switches from satu¬ 
rated to unsaturated regime. The photon number accu¬ 
mulated in the spectral peak and its observed position 
Apk is determined near Overall, we find that the 
dependence of Apk on parameters is rather weak and no 
fine-tuning is necessary to bring the peak into the ob¬ 
served range around 1 MeV. 

Note that in all models presented in this paper we 
chose the conservative assumption that the jet is initially 
photon starved — we chose a low initial photon num¬ 
ber, which corresponds to a high initial Upk ~ 10 MeV. 
Then the exact initial condition is not important — it is 
quickly “forgotten” as many more photons are generated 
at larger radii. This assumption becomes inconsistent 
for relatively slow jets, which have a large Planck radius 
(B13, V13) or for fast jets with weak dissipation below 
the Wien radius. In these cases strong thermal radiation 
should be assumed at the inner boundary of the simula¬ 
tion. 

The most efficient photon production mechanism, 
which controls the observed Upk in our models, is syn¬ 
chrotron emission in the Wien zone r < i?Wj well inside 
the photopsheric radius i?*. This mechanism works when 
a non-negligible fraction of the dissipated energy is chan¬ 
neled to nonthermal leptons. The constraints on jet mag¬ 
netization in the photon produ ction zone are less restric¬ 
tive than suggested in V13 and iThompson fc Gilll (j2014ll . 
We find that a moderate magnetization £b > 10“^ is suf¬ 
ficient to generate the photon number observed in a typ¬ 
ical GRB. VI3 required Eb ~ 1 because they neglected 
pair cascades, which increase the number of synchrotron 
emitters, partially offsetting the reduction in synchrotron 
emissivity at small ep. Also, the extended heating range 
considered in this work allows more time to accumulate 
photons in the Wien peak. 

5.2. Low energy slope a 

The typical photo n index of observed GRB spectra at 
E < Upk is a ~ — 1 (|Kaneko et al.ll200^ . In the absence 
of synchrotron photon produc tion, photospheric emission 
has a much harder spectrum (jPe ’er et al.ll2006l : iGTanniosI 
I200(H BIO; Vll), unless the jet has a very small collima- 
tion angle, comparable to T”^ (iLundman et al.ll2013ll . 
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Our transfer simulations of magnetized jets with eb > 
10“^ naturally explain the observed a ^ — 1 and its vari¬ 
ations (Figure [5]). The softening of the low-energy slope 
begins near the Wien radius where the generated low- 
energy synchrotron photons are no longer Comptonized 
to the spectral peak, ending up at intermediate ener¬ 
gies. For a broad range of parameters the resulting low- 
energy spectrum is roughly flat in number of photons 
per \tiE (i.e. a = — 1 j. Similar low-energy slopes were 
recently obtained bv lThomoson fc GiU (j20l4 l in magnet¬ 
ically dominated jets. Despite the different regime, the 
spectrum formation below the peak is similar: unsatu¬ 
rated Comptonization of a low-energy photon source. 

5.3. Baryonic vs. Poynting-dominated jets 

In this work we have only considered moderately mag¬ 
netized jets. Very weak magnetization eb < 10“^ is dis¬ 
favored based on the need for efficient photon generation 
and the observed spectral shape (in particular below the 
peak). On the other hand, our results do not exclude the 
possibility of a Poynting-dominated jet at small radii, as 
the initially dominant magnetic energy could have been 
dissipated at r <C i?*. 

Radiation from magneti callv dominated jet s wa s 
recently conside re d by iThomps on fc Gilil (20W), 
iGill fc Thompson! (|2014f) and iBeeue fc Pe’CT ( 2015 1. 
Gill & Thompson envision a two-stage evolution of the 
jet beginning from a baryon-free Poynting flux followed 
by baryon loading and photospheric emission. In their 
picture, the opacity is due to electron-positron pairs 
generated by dissipation, which takes place in two 
separate episodes. The pairs also generate photons 
through cyclo-synchrotron emission. The evolution of 
radiative processes in their scenario resembles that in 
our jet models: the generation of synchrotron photons is 
followed by Gomptonization into a Band-like spectrum. 


5.4. Models for individual bursts 


We applied our radiative transfer model to three well- 
studied bright bursts, GRB 990123, GRB 090902B, and 
GRB 130427A, which show different prompt spectra. 
Successful fits have been found in all three cases, giv¬ 
ing estimates for the main parameters L, r;, eb, F(Rc)j 
E th, and Enth (Table 1). In particular, we find the jet 
magnetizations Eb = 0.01 — 0.05, and the Lorentz factors 
between 400 and 1200. The average photospheric radii 
i?* in the three bursts vary around 10^^ cm. 

These results suggest a systematic method for esti¬ 
mating the jet Lorentz factor in a larger sample of 
bursts, which is independent of another new method 
based on the reconstruction of the GeV-(-optical flash 
produced by the external blast wave at much larger 
radii ^ 10^^ cm (iBeloborodov et al.l 120141 : iVurm et al.l 
[Ml iHascoet et al.l 1201511 . For GRB 130427A we find 


Ff = 400, which within uncertainties is consistent with 
the value Fej = 350 obtained from the GeV-(-optical flash 
reconstruction (|Vurm et al.ll20ll l. For GRB 090902B we 
find Ff Ri 1200, whic h is higher than Fpj = 600 —900 used 
to fit the GeV flash ([Hascoet et al.ll20lll . Note however 
that the reverse shock in GRB 090902B is relativistic; in 
this case there is a significant uncertainty in the upper 
limit on Fgj an d the flash modeling only gives a lower 
limit Fej > 600 ([Hascoet et al.llM(5 1. 


In GRB 130427A we find that the nonthermal heat¬ 
ing becomes weak well before the jet expands to trans¬ 
parency. Its nonthermal Band-like spectrum is mainly 
the result of thermal Comptonization. The weak resid¬ 
ual high-energy emission from nonthermal heating still 
makes a significant contribution to the GeV luminosity 
and is consistent with the Fermi LAT data during the 
main prompt emission episode. It could also explain the 
variability superimposed on the smooth extended GeV 
flash in GRB 130427A. 

The featureless Band spectrum of GRB 990123 sug¬ 
gests that thermal heating dominates the dissipation also 
in this burst. Our transfer model predicts a moderate ex¬ 
cess below a few tens of keV due to synchrotron emission, 
similar but weaker than that observed in GRB 090902B. 
A h int rf^such excesse s indeed seen in the data (Figure 

In contrast, nonthermal dissipation in GRB 090902B 
is strong up to the photosphere and beyond. It well ex¬ 
plains the observed high-energy component and the soft 
excess, which were previously model ed as a power law 
of unknown origin (jAbdo et al.ll2009H . The high-energy 
component is also expected to make a contribution to 
the GeV flash observed in GRB 090902B by Fermi LAT. 
Our result for t he Lorentz factor F 1200 is similar to 
the estimate by iPe’er et ^ (l2012fl . They used a differ¬ 
ent phenomenological model for GRB 090902B (a multi¬ 
color blackbody for the photosphere and nonthermal ra¬ 
diation from dissipation at a large radius) and estimated 
F ft! 1000. 

Confining most of the dissipation to the subphoto- 
spheric region in bursts like GRB 130427A is expected if 
dissipation is caused by neutron-proton collisions, whose 
rate declines at Tp < 20 (BIO). The heating of thermal 
electrons by Coulomb collisions with hot protons (stirred 
by n-p collisions) is also reduced at t± <C m n/(mp In A), 
wher e InA ^ 20 is the Coulomb logarithm (|Rossi et al.l 
l2006t BIO). Significant heating extending beyond the 
photosphere in GRB 090902B suggests the presence of 
a different dissipation mechanism. For instance, inter¬ 
nal shocks can occur both below and above the photo¬ 
sphere, and the shocks can directly heat the photons and 
plasma without relying on n-p or Coulomb collisions 
(IBeloborodov! [M^ . 

Note that a systematic study of the entire parameter 
space was not attempted in this paper and is a currently 
ongoing work. Thus some degeneracies may be present 
between the parameters reported in Table [U particularly 
if the actual dissipation profile is more complex than the 
“minimal” model used in this paper. The relatively fea¬ 
tureless spectra of GRB 990123 and GRB 130427A are 
more prone to such degeneracies. In contrast, the promi¬ 
nent extra component (s) in GRB 090902B make this case 
more restrictive, thus the obtained solution and jet pa¬ 
rameters are most likely unique. 

5.5. Future prospects 

Future analysis of GRB spectra using transfer simu¬ 
lations can be developed in two ways. (1) The trans¬ 
fer models give spectra that could be observed with a 
high temporal resolution dt* ^ i?*/r^c, comparable to 
a few ms. The evolving nonthermal spectrum could be 
fitted by the model; it would show the evolution of the 
jet parameters during the burst. (2) Photon statistics 
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in observed GRB spectra are usually accumulated on 
timescales 5t ^ i 5G, hindering the resolution of the in¬ 
stantaneous spectrum emitted by the variable jet. Then 
even the best time-resolved data analysis may give a mix¬ 
ture of different instantaneous spectra, which could result 
in the presence of multiple components in the measured 
spectrum. Transfer simulations may be used for the anal¬ 
ysis of the multiple components. A dominant Band com¬ 
ponent is almost always found in GRBs, including the 
recent Hts by the evolving mix ture of three compo nents: 
Band -I- power law -|- thermal (|Guiriec et al.l[2(Tl5l l. The 
nonthermal emission may be well explained by our model 
of a heated jet, including the excess at low and high en¬ 
ergies which was previously viewed as a separate power 
law (see Figure [T]). The model also predicts that emis¬ 
sion from weakly heated jets has a quasi-thermal shape. 
Thus the presence of a quasi-thermal component (iRvdel 


12004 [2005t iRvde fc Pe’erl[2009l : iGuiriec et ani2015ll may 

indicate the presence of unresolved parts of the jet with 
weak heating. 

The transfer model makes specific predictions for po¬ 
larization of the prompt radiation (Lundman et al. 2016, 
in preparation), which can be tested by future ob¬ 
servations. The polarized radiation arises from syn¬ 
chrotron emission by nonthermal electrons; it increases 
toward lower frequencies below the MeV peak and is 
strongest in the X-ray band. Bursts with strong non¬ 
thermal dissipation extending to the photosphere (such 
as GRB 090902B) will be most promising for the detec¬ 
tion of polarization. 

We thank Ghristoffer Lundman for helpful comments 
and discussions that helped to improve this manuscript. 
This work was supported by NSF grant AST-1412485 
and NASA ATP grant NNX15AE26G. 


APPENDIX 

A. DISSIPATIVE ,IET DYNAMICS 


Let us rewrite the radiative transfer equation (O in the form 

r^F^ 51nr ^ ~ ~ + (1 + m)(1 ~ d) 

f) T f) 

+ (1 + M)(1 - ~ 

Integrating Equation over the photon energy and taking the first two angular moments yields 



(A2) 


(A3) 


where j is the angle-averaged bolometric emissivity and we have assumed an isotropic (comoving) opacity k. The 
moments of intensity Ii, I 2 and I 3 are defined as 


-fm = ^ J di^ J d^i. 


(A4) 


The lab frame luminosity is given by Lrad = (4 7r)^r ^F^(/ o -|- 2 /i + 12 ) (using F ^ 1, see Bll). The equation governing 
its evolution is obtained by adding Equations (IA2I) and (IA3|) . 


dLrad 

dlnr 


{4TTfr^r[j-K{Io + Ii)] . 


(A5) 


The dissipation rate in the comoving frame is dE\^/[dV'dt') = 47r(j — kIq). Transformation to the lab frame yields 


Equation (IA5I) thus becomes 


dlnr dV'dt' 


{4-Kfr'^T (j - kIq). 


(A6) 


dLj-ad 
din r 


dlnr 


(47r)^r^F kIi. 


(A7) 


The second term on the right hand side describes the work done by the radiation field on accelerating the jet. 

If dissipation is described as redistribution of the total energy L = Lp\ + Trad = const from the bulk kinetic (Lpi) 
to internal (radiation-dominated) form, one can write Trad = L(1 — F/yy) to obtain a dynamical equation for the jet 
Lorentz factor 


dr 


de-k 

= -V-r- + (^tZ± 
dr 


AttIi 

rripC^ 


(AS) 
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where dsh/dlnr = L~^dL\^/dh\r. Alternatively, the heat source can be modelled as a “reservoir” of free energy carried 
by the jet in the form of internal bulk motions, iturb- Then L = Lp\ + Lrad + Tturb = const, and the energy balance 
is given by 


dL\i _ dlj^ad ^ 

d In r d In r d In r 


(A9) 


In this case, the dynamical equation becomes 


dr 




47r/i 

rripC^ 


(AlO) 


In the optically thick domain it is useful to rewrite the dynamical equation in a form where the (small) first moment 
Ii o f the intensity does not appear. In such case I 2 ~ /q/3 » Ii and L^ad = ( 4 / 3 )( 47 r)^r^r^/Q, whereby Equation 
(|A2I) can be written as 


dLrad 

dlnr 


4 dAh 
3 dlnr 



where the last term accounts for adiabatic cooling. In place of Equation (IA 8 I) one now obtains 


(All) 


dluE 1 
dlnr 2 + r /?7 

In the absence of dissipation Equation (IA12|) can be int egrated straightfo rwardly to yield the standard solution for E 
in baryonic (initially) radiation dominated jets (see e. g.lPiran et al.iri993fl . 

In the reservoir model of the heat source Equation (lAlOl) is replaced by 


2 ( 1 - 1 ] 

4 

. V vj 

dlnr 


(A12) 


dlnE 2 Arad ~ dLh/dlnr 
dlnr 3TMc^ + 2Lrad 


(A13) 


B. PAIR-LOADING IN A DISSIPATIVE JET 

The cooling of the injected pairs and their subsequent thermalization is very efficient in the high compactness 
environment of GRB jets. Thus bulk of the pairs reside in a Maxwellian distribution at any given time. Their 
density is controlled by injection of new pairs (both primary particles as well as secondaries from the cascade), 
annihilation, and expansion of the jet (cf. BIO), 


where 


c d 

r2 fij. 


(rr^n±) 


— Tl^ ^anri; 


fl± 


dTTmeC^r^r ’ 


^ann 


3 3 

-cc7Tn+n_ = —cdT 

4 16 


(nl-nl), 


(Bl) 


(B2) 


and n± = + n_ = 2n+ + rip. The pair yield Y characterizes the fraction of energy injected into the primary pairs 

that ends up in pair rest mass through the pair-photon cascade. Defining the pair loading Z± = n±jnp, and using 
Tr^rip = M/( 47 rmpc) = constant. Equation (IBII) can be written as 


dZ^ _ YcuthL TUp 3 172 -1 \ 

dlnr “ rMc2 rUe 


(B3) 


where Tp = CTripr/r. _ 

Two regimes can be identified in Equation (IB3I) : (1) creation-annihilation balance, where both terms on the RHS are 
larger than the LHS, and (2) freezeout, where annihilation (last term) is negligible. In creation-annihilation balance 
the pair loading evolves as 




f IbEenthA Trip 
ySrpTMc^ me 



(B4) 


Annihilation freezes out once t± = Z±Tp = 16/3; beyond this radius the pair loading is governed by the first two terms 
in Equation (IB3I) and can increase at most by a logarithmic factor in radius (assuming a logarithmically flat heating 
rate and F = const). 


C. PHOTON PRODUCTION BY SYNCHROTRON EMISSION 

The rate at which photons are generated and accumulated in the spectral peak depends on the interplay between 
synchrotron emission, reabsorption, IC scattering and induced Compton downscattering. It was shown in VI3 that 
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competition between these processes defines a characteristic photon energy above which the emitted synchrotron 
photons are upscattered to the Wien peak rather than downscattered and/or reabsorbed. Below we extend the results 
of V13 to include pair-photon cascades and apply them to the extended heated Wien zone. 

At high optical depths the timescales for all relevant radiative processes are much shorter than the jet expansion 
timescale. In this regime the Kompaneets equation can be written as 


CTTne 1 _ ^^4 

ITleC dE 


oE 


= jv{E) — CK{E)n, 


(Cl) 


where n is the photon occupation number, is the emissivity (here we consider only synchrotron) in units of the 
number of photons per unit phase space volume per second, k is the (synchrotron absorption) opacity, rie is the density 
of the thermal electron/positron component, and is the electron temperature. We have neglected the recoil term in 
Equation (ICll) since photon generation t akes place at A <C k-^Te <C me(? where recoil losses are negligible. 

By comparing the terms in Equation (ICll) one can identify three regimes where different physical processes domi¬ 
nate. At low enough photon energies the radiation has the usual optically thick synchrotron spectrum for power-law 
electrons 0 


CK p+2 \Eb J \meC^ ) 


(C2) 


As the synchrotron emitting electrons are relativistic, the brightness temperature of the optically thick synchrotron 
radiation Tb = En/kB ^ Tg. In this case the induced scattering term dominates the Kompaneets operator, tending to 
downscatter the synchrotron photons. The first break in the spectrum occurs where the induced scattering rate becomes 
dominant over the reabsorption rate. Above this energy (but before the second break, see below) the Kompaneets 
equation reduces to 


<JTng 1 a 4 2 

rUeC E'^ dE 


ME). 


(C3) 


The synchrotron emissivity of a power-law distribution of electrons, n.e(7) = n-o"/ is 


>(£’) = jo 


/ \ -(P+5)/2 


^ -1 

—ttf caTno 


f rrigC^ 

V 


2 


f E \ 


(C4) 


where af 


jhc if the fine structure constant and Eb = 


h eB/nigC. The solution of Equation (IC3I) takes the form 


n = 


2meCj„{E) 


1/2 


{p - l)aTneE 


^ E-iP+^)/\ 


(C5) 


Compton upscattering begins where the radiation brightness temperature decreases down to the electron temper¬ 
ature. This yields the second break in the spectrum, above which the usual thermal Comptonization spectrum is 
established. The characteristic energy, Eq, can be found by equating the two terms in the Comptonization operator 
in Equation (ICll) . 


Using the spectrum (IC5I) . this yields 


9ri _ _2 

Me ■ 


or equivalently 


p + 7 kBTg 


n{Eo), 


(jp- l)(p-|-7)^ (/cBTe)^ gT^e 
32 Uo rrigC 


ji^iEo). 


(C6) 


(C7) 


(C8) 


To proceed, one has to specify the electron distribution. Let us assume that relativistic electrons are continuously 
injected at a fixed Lorentz factor yinj with power Qinj [erg • cm“^ • s“^]. The injected high-energy pairs initiate a 
pair-photon cascade with a large number of secondaries. Let us denote their number per injected primary particle as 
Mg. Denoting the pair yield as K = Als/7inj we can express the equilibrium pair distribution in the cascade as 


Ml) = nol 


(C9) 


^ Here we are using the delta-function approximation, assuming 
that all synchrotron photons from a single electron are emitted at 


E = 0.3^'^Eb- 
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where 


no = 


^Yl±. 


min 


■S^nth 


cT 


dcCTT (SB + Erad) r 


(CIO) 


p = 2 + a±, 7 ±,min is the electron energy below which the saturated cascade turns off, and a± ~ In A^s/(ln 7 inj — 
ln 7 ±,min)- The energy fractions £b and £rad are dehned via 


17b = 


£bL 


and Urad = 


EradT 


(Cll) 


47rcr2 dTTcr^r^ ’ 

and £nth is the fraction of total available energy deposited into the injected electrons per logarithmic radius interval, 


EnthT = 

alnr 


mj- 


Using riQ given by Equation (IC10|) in the expression (IC4I) for the emissivity and inserting the latter into 
the critical photon energy/electron Lorentz factor 


Eo \ 


= (703 70 )^+' = 


967r 


EnthU7^ 

,min tt f Eb 


(p-l)(p + 7)2 ^ (EB+Erad) y'^\meC'^ 


Here tt = aTZ±npr/T, y = is the Compton parameter and the baryon number density is given by 

L 


AirnipC^r^Tr] 


(C12) 

yields 

(C13) 

(CM) 


The value of the critical Lorentz factor 70 is typically between a few and a few tens. If the pair cascades are suppressed 
(e.g. by synchrotron cooling), one obtains 


^ - 7 QT -3/10 n2/5 p3/5 1/5 -1/5 ^-2/5 -1/10 

70 — 7.0^52 Tiio 1 % 7inj.2 ^-2 ^B 

The synchrotron photon production rate at E > Eq is 

.^Igynch — 


Enth 


Eb + Erad 


1/5 


47r 

{chy 


f liE)E^dE= ^ gnthy7±fndn Cr^ 

Eg 2(p — 1) (EB+Erad) T Eb \Eb, 


At a sphere of radius r, the number flux of photons produced between rmin and r is 

tVph)?') = [ 47rr^A4ynch(£)dlnr. 


(C15) 


(C16) 


(C17) 


The quan tities Erad(i') and r(r) are determined by the heating model, in particular by dL^/dlnr, see Equations 
(|A11I) and (IA13I) . One last quantity entering the photon production rate is the electron temperature. In the Wien 
zone it can be determined from the relation 


Trad — ^EradT — dr/uB+e Eph- 


(CIS) 


For given cascade parameters Y and 7 ±,min, Equations (|CI3|) - (ICI 8 I) along with (lAIIl) and (IAI3I) form a closed set. 
Its solution gives the number of synchrotron photons 77ph(£) which are Comptonized toward the Wien peak. The 
solution is simplified if one uses the approximation Np^ ~ A4ynch • 

It is convenient to express the result as the number of produced photons per proton, riph/up = Npi^/Np, where 
Np = L/^mpC^r]) is the proton number flux carried by the flow. In the absence of cascade, q!± = 0 and Y = 7 “^ 


one hnds 


Uph 


= 4.1 X 10'’T 


6r-i/7 3/7p-5/7 6/7 -4/7 3/7 2/7 


52 


' 10 


% 7inj,2 ^B ^rad 


^nth 


H” ^rad 


4/7 


In the opposite limit of a fully saturated cascade one obtains 


^ = 1.83 X 10^ rr^ sill I 


± Vs ^B ^rad 


Eb + Erad 


2/5 


(CI9) 


(C20) 
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where we used p = 3 {a± 


1). In a general case, the radial dependence of photon production is given by 


^ph 


3p-|-l p 3p-l-l 


p— 1 2p —2 

735+T _ 35 TT 
^rad 


/fnthr7^y'’+ 

y £B + £rad j 


(C21) 


The above analysis neglects some effects that are not straightforward to describe analytically. In particular, near 
the critical Lorentz factor 70 ^ 10, the IC cooling is modified by the Klein-Nishina effect, as the target photon energy 
in the electron rest frame ^ SfceTeyo is comparable to m^c^. The Klein-Nishina recoil effect also influences the fate 
of the IC photons — they can scatter and lose energy to electron recoil before annihilating into pairs. This effect 
is particularly important for those IC photons that see a low opacity for pair production, which happens if they are 
below the threshold for interacting with the Wien-peak radiation. These effects can substantially alter the number of 
electrons near 70 , and consequently the number of produced synchrotron photons. The above analytic estimates can 
only serve as a rough guide to the photon production showing its trends, i.e. the dependence on r, £b, F etc. The 
accurate photon number is provided by our numerical simulations. 
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